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> 

■ Bose-Einstein correlations of charged and neutral kaons have been measured in 

^ ■ e ± p deep inelastic scattering with an integrated luminosity of 121 pb _1 using 

the ZEUS detector at HERA. The two-particle correlation function was stu- 
died as a function of the four-momentum difference of the kaon pairs, Q12 = 
y— (Pi — P2) 2 5 assuming a Gaussian shape for the particle source. The values of 
the radius of the production volume, r, and of the correlation strength, A, were 
obtained for both neutral and charged kaons. The radii for charged and neutral 
kaons are similar and are consistent with those obtained at LEP. 
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1 Introduction 



The use of Bose-Einstein correlations (BEC) in particle physics as a method of deter- 
mininig the size and the shape of the source from which particles originate was first 
considered by Goldhaber et al. [1, 2] in 1959 for pp annihilations. Bose-Einstein corre- 
lations originate from the symmetrization of the two-particle wave function of identical 
bosons and lead to an enhancement of boson pairs emitted with small relative momenta. 
The effect is sensitive to the size of the emitting source. The studies of BEC for pairs 
of identical particles have been carried out in a large variety of particle interactions. In 
particular, HI and ZEUS have reported results on inclusive charged particle pairs in e p 
collisions at HERA [3,4]. Recent reviews [5-9] summarise the underlying theoretical con- 
cepts and experimental results. The measurements of the radius of the emission source 
have been mostly performed for neutral and charged pions. For other bosons, e.g. kaons, 
the information is scanty. 

This paper reports first results on BEC for pairs of charged (K ± K ± ) and neutral (K^Kg) 
kaons in deep inelastic scattering (DIS) at HERA. The measurements are compared with 
results from e + e~ interactions whose fragmentation properties are expected to be similar 
to the current region of DIS [10]. However, proton fragmentation may lead to a significant 
difference in the properties of the hadronic final state. 

The correlation function for two identical kaons is defined as 



where Qu is the four-momenta difference of the kaons with four momenta p l and p 2 given 

as 



where Mkk is the invariant mass of the pair of kaons and rriK is the kaon rest mass. The 
function P{Qi2) in Eq. (1) is the two-particle density: P{Q\2) = (^/N)(driKK /dQn) , 
where ukk is the number of kaon pairs and N is the number of events. The denominator 
of Eq. (1), Po^Qu), is the two-particle density in the absence of BEC. 

For a static source with a Gaussian density distribution, the correlation function can be 
parametrised as follows [2]: 



P(Ql2) 



(1) 




(2) 



R{Q 12 ) = l + \exp(-r 2 Q 2 12 ). 



(3) 
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The A parameter gives information about the strength of the BEC. For a completely 
coherent source, A is zero, while for a completely incoherent source, A is one [11]. Con- 
tributions from decays of short-lived resonances can further modify this parameter. The 
parameter r is related to the size of the source and is called the radius in the following 
sections. 

2 Experimental set-up 

The analysis was performed with data taken by the ZEUS detector between 1996 and 2000 
at HERA. The data from e^p collisions collected in this period with electron 1 energy E e = 
27.5 GeV and proton energy E p = 820 GeV (1996 - 1997) or E p = 920 GeV (1998 - 2000) 
correspond to an integrated luminosity of 121 pb" 1 . 

A detailed description of the ZEUS detector can be found elsewhere [12]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [13], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD con- 
sists of 72 cylindrical drift chamber layers, organized in 9 superlayers covering the polar- 
angle 2 region 15° < 9 < 164°. The transverse-momentum resolution for full-length tracks 
is o~(pt)/pt = 0.0058pT © 0.0065 © 0.0014/pt, with px in GeV. To estimate the energy 
loss, dE/dx, of tracks, the truncated mean of the sense-wire pulse-heights was recorded 
for each track, discarding the 10% lowest and up to the 30% highest pulses [14-16]. The 
measured dE/dx values were normalised to the dE/dx peak position for tracks with mo- 
menta 0.3 < p < 0.4 GeV, the region of minimum ionisation for pions. Henceforth, 
dE/dx is quoted in units of minimum ionising particles (mips). The resolution of the 
dE/dx measurement for full-length tracks is about 9%. The tracking system was used to 
establish the primary and secondary vertices. 

The high-resolution uranium-scintillator calorimeter (CAL) [17] consists of three parts: 
the forward, the barrel and the rear (RCAL) calorimeters. Each part is subdivided trans- 
versely into towers and longitudinally into one electromagnetic section and either one or 
two hadronic sections. The smallest subdivision of the calorimeter is called a cell. The 
CAL energy resolutions, as measured under test-beam conditions, are <r(E) / E = 0.18/ \f~E 
for electrons and a(E)/E = 0.35/^/E for hadrons, with E in GeV. 



1 Here and in the following, the term electron" denotes generically both the electron and the positron. 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the center of HERA. The coordinate origin is at the nominal interaction point. 



2 



The energy of the scattered electron was corrected for energy loss in the material between 
the interaction point and the calorimeter using the small-angle rear tracking detector 
[18,19] and the presampler [18,20]. 

3 Event and track selection 

The inclusive neutral current DIS process e(k) + p(P) — > e(k') + X can be described in 
terms of the following kinematic variables: Q 2 , the virtuality of the exchanged photon, 
x, the Bjorken scaling variable and y, the fraction of the lepton energy transferred to the 
proton in the proton rest frame. They are defined as follows: Q 2 = —q 2 = — (k — k') 2 ; 
x = Q 2 /(2P ■ q); y = (q ■ P)/(k ■ P), where k, k' and P are the four-momenta of initial 
and final scattered electrons and incoming proton, respectively. These variables were 
reconstructed using the electron method (denoted by the subscript e), which requires 
measurements of the energy and angle of the scattered electron. 

A three-level trigger system [12] was used to select events online. At the third level, 
electrons with energy greater than 4 GeV and position outside a rectangle defined by 
|X| < 12 cm, |Y| < 6 cm on the face of the RCAL were accepted. Data below Q 2 ~ 
20 GeV 2 were prescaled to reduce the trigger rate. 

The offline selection of DIS events was based on the following requirements: 

• |Z vtx | < 50 cm, where Z vtx is the Z component of the primary- vertex position deter- 
mined from the tracks. This cut reduces the background from non-ep interactions; 

• an identified scattered electron in the CAL with energy E e > 8.5 GeV; 

• 2 < Ql < 15000 GeV 2 ; 

• 35 < 5 < 60 GeV, where 5 = ^2^(1 — cos 9i), E^ is the energy of the i th calorimeter 
cell and 9i is its polar angle as viewed from the primary vertex. The sum runs over all 
CAL cells. This cut reduces the background from photoproduction and events with 
large initial-state radiation; 

• ye < 0.95, to remove events with misidentified scattered electrons; 

• 2/jb > 0.04, to remove events with low hadronic activity, where i/jb is the value of y 
reconstructed using the Jacquet-Blondel method [21]. 

Good quality tracks measured in the CTD with high acceptance and resolution were 
selected using the following requirements: transverse momentum px > 0.15 GeV and 
pseudorapidity \r)\ < 1.75. In addition, the tracks were required to pass through more 
than three CTD superlayers. 
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After the above cuts, the data sample contained 25 million events with at least two good 
tracks and the average Q 2 of the sample was (Q 2 ) = 35 GeV 2 . 

3.1 Selection of charged kaons 

Charged kaons were selected using the energy-loss measurement, dE/dx. The analysis 
used all tracks fitted to the primary vertex with the exception of the scattered-electron 
track. Tracks were selected as described above. The dE/dx as a function of momentum 
for positively charged tracks is shown in Fig. la. The curves indicate the region used for 
identification of positively charged kaons. Kaons were selected by requiring / < dE/dx < 
E, where / and F are functions of the track momentum, p, motivated by the Bethe-Bloch 
equation. For positive kaons / = 0.08/p 2 + 1.0, F = 0.17/p 2 + 1.03 mips and for negative 
kaons / = 0.08/p 2 + 1.0 and F = 0.18/ p 2 + 1.03 mips (with p in GeV). The slight 
difference arises from the different response of the CTD to positive and negative tracks. 
The kaons were identified for p < 0.9 GeV and dE/dx > 1.25 mips. 

The kaon identification efficiency for p t > 0.15 GeV, \rj\ < 1.5 and p < 0.9 GeV was 61%, 
with a purity of 90%. The resulting data sample contained 55522 K + K + or K~K~ pairs. 

3.2 Neutral kaon selection 

The K s mesons were identified using the charged-decay channel K s — > 7r + 7r~ with a 
similar selection as in a previous publication [22]. The pion tracks were required to 
originate from secondary vertices. Assigning the pion mass to both tracks, the invariant 
mass M(n + n~) was calculated and the candidate was accepted if the mass was within 
± 20 MeV of the nominal PDG [23] K s mass. To eliminate tracks from photon conversions 
and A/A contamination, the electron, pion and proton masses were assigned to tracks and 
the following cuts were used: M(e + e~) > 80 MeV and M(np) > 1121 MeV. 

The following additional requirements were applied to the selected K s candidates: 

• 2 < Ld < 30 cm, where Ld is the decay length of the K s candidate; 

• AZ < 0.8 cm, where AZ is the projection on the Z axis of the vector defined by the 
primary interaction point and the point of closest approach of the K s candidate; 

• «xy < 8°, where «xy is the (collinearity) angle between the candidate Kg momentum 
vector and the vector defined by the interaction point and the K s decay vertex in the 
XY plane; 

• pf A > 0.11 GeV, where the Podolanski-Armenteros variable pf A is the projection of 
the candidate pion momentum onto a plane perpendicular to the K s momentum di- 
rection [24]. 
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The total number of Kg candidates was 725505. After all cuts, the selected data sample 
contained 19494 K S K S pairs and 400 triples. Each combination of two particles was 
included in the analysis. Figure lb shows the 7T + ir~ invariant mass distribution after the 
J^-pair selection and a fit to the signal plus linear background, which resulted in an 
estimated background of 1.4%. 

4 Monte Carlo simulation 

Inclusive DIS events with Q 2 > 2 GeV 2 were generated without BEC using the Ariadne 
4.10 Monte Carlo (MC) model [25] interfaced with Heracles 4.6.1 [26] via the Djangoh 
1.1 program [27,28] in order to incorporate first-order electroweak corrections. The Lund 
string model [29] was used for the description of hadronisation, as implemented in the 
Jetset 7.4 [30] program. 

The generated events were passed through a full simulation of the detector using the 
Geant 3.13 program [31] and reconstructed and analyzed in the same way as the data. 
The Ariadne MC sample corresponds to a similar integrated luminosity as that of the 
data. 

5 Extraction of BEC parameters 

The main difficulty in measuring BEC is in the construction of a reference sample which 
should be as close as possible to the analyzed data in all aspects but free from the Bose- 
Einstein effect. The obvious reference sample provided by unlike-sign charged kaon pairs 
cannot be used due to the strong signal of the 0°(1O2O) — > K + K~ decay at low values of 

Ql2- 

A reference sample can be derived from a Monte Carlo simulation without BEC. In this 
so-called single-ratio method, the correlation function R(Q 12) is defined as: R s (Qu) = 
-P(Qi2) data /-P(<5i2) MC ' noBEC , where P(Qi2) data is the normalized two-particle density dis- 
tribution for the data and P(Qi2) MC ' noBEC is the corresponding distribution obtained for 
MC without BEC. However, this approach requires a correct simulation of the physics 
processes in the absence of BEC, as well as a good description of the detector effects. 

In another approach, a reference sample can be obtained using an event-mixing procedure 
where two kaons from different events are combined. This method, which is used in 
this analysis, is less sensitive to imperfections in the MC simulation. To correct for 
other correlations lost in the event-mixing procedure, the two-particle correlation function 
R{Qi2) w &s calculated using the double-ratio method: 



5 



R(n , P(Qi 2 ) data I P(Q 12 ) MC '™ BEC m 

Wl2) p mix (gi 2 ) dat 7 p m ix(Qi 2 ) MC ' noBEC ' 1 ] 

where -P m ix(Qi2) data is the two-particle density constructed from pairs of kaons coming 
from different events and -P m ix(Qi2) M oBEC is obtained in a similar way for MC events. 
The double-ratio method was used for the main analysis and the single-ratio method only 
to estimate systematic uncertainties. 

To fit the correlation function R(Qu) defined by Eq. (4), the modified Goldhaber parametri- 
sation (Eq. (3)) multiplied by an empirical term 1 + /3Q 12 , which accounts for the presence 
of possible long-range two-particle correlations for high Q12, is often used: 



R(Qi2) = «(1 + Ae-^ r2 )(l + (3Q l2 ). (5) 

Such correlations are imposed for example by energy and charge conservation, phase-space 
constraints or strangeness compensation. In this analysis, the f3 parameter was found to 
be zero within errors and its possible deviation from zero was included in the systematic 
uncertainties. 



6 Systematic uncertainties 



Systematic uncertainties on the fitted A and r parameters arise from event and track 
selection, the modeling of dE/dx, the fitting procedure and the construction of the ref- 
erence sample. They were calculated from the deviation of the fit parameters from their 
nominal values after changing the analysis cuts or procedures. For charged kaons, a bias 
originating from the contamination of the experimental kaon sample by pion, proton and 
antiproton is expected. For neutral kaons, a bias can be introduced by the contamination 
of the sample by Kg Kg pairs from the decay of /o(980) resonance. These effects are 
discussed in the next sections. 

The following systematic studies were carried out for the charged kaon sample. The 
resulting changes for A and r are given in parentheses as [AA, Ar]: 

• the correlation function was calculated using the single-ratio method [—0.01, —0.03]; 

• the fit was repeated for different lower and upper limits of Qi2- In addition the binning 
in Q 12 was modified [l^t-ao' ]; 

• the momentum range was reduced to p < 0.7 GeV for both the double [+0.06, +0.06] 
and single [+0.01, —0.01] ratio methods; 
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• the definition of the kaon band was varied for the double-ratio method [—0.02, +0.05] 
and the single-ratio method [—0.04, —0.02]; 

• the track quality cuts were changed within their resolutions [—0.02, +0.04]; 

• the DIS selection cuts on E e , y e , t/jb, 5 were varied within the resolutions [lo!o2)^o!o2 ]> 

• the influence of pion, proton and antiproton contamination was checked by tightening 
and relaxing the kaon selection criteria. The purity was raised to 99% [+0.05, +0.07] 
and was brought down to 80% [—0.04, —0.05]. 

The following systematic uncertainties for neutral kaons were considered: 

• the correlation function was calculated using the single-ratio method [+0.10, +0.02]; 

• the fit was repeated for different lower and upper limits of Qi2- In addition, the binning 
in Q 12 was modified [1^+^]; 

• the cuts on the K s momentum, axY> AZ, M(ir + iT~) were varied [—0.08, —0.08]; 

• the cut on pf A was changed to 0.12 GeV [+0.04, +0.02]; 

• the full Eq. (5) was used to check the sensitivity of the fit to possible long-range 
correlations [+0.11, +0.03]; 

• the track quality cuts were changed within the resolution [+0.15, +0.02]; 

• the cut to remove the e + e~ background was changed to M(e + e~) > 50 MeV [+0.08, +0.03]; 

• the DIS selection cuts on E e , y e , yjB, 5 were varied within resolutions [laoL-aoiH- 

The contributions from the different groups of systematic uncertainties in the parameters 
A and r were added in quadrature separately for positive and negative variations. The 
overall systematic uncertainties in A and r for K ± K ± and K S K S are presented in Table 1. 



7 Results 

7.1 Correlations in K ± K ± pairs 

Figure 2 shows the measured two-particle correlation function for K K ± pairs. The 
results obtained by the fit function given by Eq. (5) with (3 = are r = 0.57 ± 0.09 
fm and A = 0.31 ± 0.06. The 90% purity of the kaon selection introduces in the kaon 
pair sample an 18% admixture of unlike particles pairs which are not correlated. The 
measured A is therefore expected to be underestimated by 18%. The systematic checks 
involving purity confirm this expectation. After the correction for purity, the result is 
A = 0.37 ± 0.07±g;g|. The value of the radius is not affected by this correction and 
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the result is r = 0.57 ± 0.09^Q g fm. The corrected parameters of BEC correlations 
for K ± K ± pairs are presented in Table 1. The radius value is consistent with that for 
charged particles: r = 0.54 ± 0.03+°;^ fm for HI [3] and r = 0.666 ± 0.0091°;°^ fm 
for ZEUS [4]. The obtained value of A for kaons is in agreement with the HI result for 
charged particles A = 0.32 ±0.02 ±0.06, although somewhat smaller than the ZEUS result 
A = 0.475 ± 0.00712 nm- 



7.2 Correlations in KgKg pairs 

The KgKg pairs may originate not only from K°K°, K°K° (strangeness S = ±2) states 
but also from the K°K° (strangeness S — 0) system, which may come from the decay 
of resonances. It has been shown [32, 33] that a Bose-Einstein-like enhancement is ne- 
vertheless expected in the Qu distribution of the KgKg pairs, even when their origin is 
from the K°K° system. According to the MC simulation of DIS events, about 75% of 
I0W-Q12 KgKg pairs come from K°K°. Similar to the case of charged kaons, Eq. (5) with 
(3 = was used to fit the correlation function calculated for K° S K^ pairs. The so-called 
raw results are A = 1.16 ± 0.29lQ;Qg and r = 0.61 ±0.081°;^ fm. The results are presented 
in Table 1 and Fig. 3. The measured radii for K ± K ± and KgK® are close to each other. 

For KgKg the fit does not take into account a possible contamination from the scalar 
/o(980) resonance decaying below the KK threshold. The decay channel /o(980) — > K°K° 
is not included in the standard MC, therefore the / (980) contribution may distort the 
strength of the Bose-Einstein effect. In order to investivate this, a MC sample with BEC 
at maximum strength was generated. A comparison of this MC with the data showed an 
excess for Q12 < 0.6 GeV, which may indicate a possible fo contribution. The same excess 
is also visible in the M^k distribution, as illustrated in Fig. 4a. 

The expected contribution from f can be approximated by a modified Breit-Wigner 
function as proposed by Flatte [34-36]: 

da N m l ■ T kk , fi x 

' 7~2 T72 vTTTZ rF 7~n v\a > W 



dM KK {rnl-Ml K y + {m Q -{T^ + r KK )Y 

where m = 0.954 GeV is the mass of /o(980) and the widths T n7T and Tkk are related 
to the coupling constants = 0.11 and gx = 0.423 by = g^ M KK / 4 — and 
Tkk = 9k^/M kk /A — m K . The normalization factor, Np, was adjusted to give the total 
number of /o(980) — > KgK® decays for Mkk < 1.8 GeV. Fig. 4a shows the distribution 
of da/dM KK according to Eq. (6). Figure 4b shows the M KK distribution in the data 
after subtraction of the standard MC without BEC together with different contributions 
of the /o resonance. An f contribution of about 7% is sufficient to account for the data, 
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therefore, the present analysis cannot distinguish the Bose-Einstein effect from the effect 
of an /o contribution. Assuming that both BEC and fo are present, the amount of fo was 
estimated from the shape of the R(Q 12) distribution (Eq. (4)). For this purpose, fits of 
the correlation function were carried out with different percentages, c/ , of fo subtracted 
from the data. A shallow minimum of the X 2 ( c /o) value for Cf — 4% was found, with 
one-sigma limits of c/ = 1% and c/ = 7%. The values of A and r at c/ = 4% were taken 
as the most probable result and their uncertainties were calculated from the one sigma 
limits of Cf . 

The results corrected for the fo contamination are included in Table 1. After adding 
the uncertainties coming from the fo subtraction to systematics, the corrected results 
are A = 0.70 ± 0. 19^53 and r = 0.63 ± O.OOl^og fm. The uncertainty on a possible fo 
admixture leads to a large systematic uncertainty on A. The radius r is less sensitive to 
the fo admixture. The radii for K^Kg and K ± K ± are very similar. 

7.3 Comparison with e + e interactions 

Figure 5 shows the comparison for r between DIS and e + e~ annihilation results at 
LEP [36-39] for both charged and neutral kaons. The radius value obtained in DIS agrees 
with the measurements from LEP. The same figure also presents DIS results for uniden- 
tified charged particles. The kaon results agree with those for charged particles within 
systematic uncertainties. The A value for charged kaons in DIS is somewhat smaller than 
that in e + e~ collisions, which may be related to the fact that the DIS data mostly populate 
the proton fragmentation region in the Breit frame. 

8 Conclusions 

Bose-Einstein correlations have been measured for pairs of charged and neutral kaons in 
deep inelastic scattering at HERA using the ZEUS detector. The values of the radius 
for charged and neutral kaons agree within systematic uncertainties. They are also con- 
sistent with the measurements for charged particles in DIS and kaons in e + e~ collisions 
at LEP. The /o(980) — > KgKg decay can significantly affect the A parameter for K^K° S 
correlations. 
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Figure 1: (a) The energy-loss dE/dx as a function of the momentum p for tracks 
with positive charge. The tracks with f < dE/dx < F, dE/dx > 1.25 mips and 
p < 0.9 GeV were taken as K + . (b) The ir + ir~ invariant-mass distribution of the 
Kg candidates. The solid line shows the result of a fit with the sum of two Gaussian 
functions and linear background. The two Gaussians were fixed to the same mean 
value. A mass in agreement with the Particle Data Group value [23] and a width 
of 3.2 MeV (6.3 MeV) were obtained from the fit for the central (second) Gaussian. 
The dashed line shows the linear background. The short vertical lines show the 
mass window used to define the signal. 
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Figure 2: T/ie two-particle correlation function for charged kaons with a fit to 
Eq. (5), with (3 = 0. The error bars represent the statistical uncertainties. 
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Figure 3: The two-particle correlation function for neutral kaons with a fit to 
Eq. (5), with f3 = 0. The error bars represent the statistical uncertainties. 
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Figure 4: (a) KgKg invariant mass for data (points) and Monte Carlo with BEC 
(dotted histogram) and without BEC (full-line histogram). The dashed line shows 
the Flatte function used to describe the fo resonance, (b ) The difference between 
data and MC without BEC compared to different fo contributions. 
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Figure 5: Comparison of DIS and LEP results for r obtained from Bose-Einstein 
correlation studies of charged and neutral kaons and unidentified charged particles. 
The DIS result for neutral kaons is corrected for the fo contribution. 



1(3 





A 


r [fm] 


K ± K ± (corrected) 

K° S K° S (raw) 
K%K% (corrected) 


0.37 ±0.07 t° fs 
1.16 ±0.29 
n 70 ±019 +°- 28 + - 38 

U.IU 3= U.1J _ Q 08-0.52 


0.57 ±0.09 
0.61 ± 0.08 t°oZ 

o.63 ± o.o9 ±8:g±8:g 



Table 1: TTie radius r and the correlation strength A /or charged and neutral kaons 
extracted from fitting the Goldhaber parametrization (Eq. (5) with (3 = 0) to the 
Bose-Einstein correlation function. The first uncertainties are statistical and the 
second systematic. The raw results for KgK® pairs correspond to those shown in 
Fig. 3; the corrected results were obtained after subtraction of the /o contribution 
from the data. For the corrected values the third uncertanty comes from the f 
subtraction. 
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